Maintaining the local temperature below the critical value in thermally out of 
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A generalized theory of open quantum systems combined with mean-field theory is used to study 
a superconducting wire in contact with thermal baths at different temperatures. It is shown that, 
depending on the temperature of the colder bath, the temperature of the hotter bath can greatly 
exceed the equilibrium critical temperature, and still the local temperature in the wire is maintained 
below the critical temperature and hence the wire remains in the superconducting state. The effects 
of contact areas and disorder are studied. Finally, an experimental setup is suggested to test our 
predictions. 
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Ever since the discovery of superconductivity, fabricat- 
ing a superconducting (SC) wire that conducts electricity 
without dissipation at room temperature has been a ma- 
jor goal of modern condensed matter physics. However, 
since the discovery of High- T c superconductors [1] (which 
have T c as high as ^150K, still far from room tempera- 
ture and perhaps close to the upper limit in these mate- 
rials [2] ) there has been little or no progress in increasing 
T c . Recently, several suggestions have been put forward 
to increase T c by fabricating nanoclusters [3, 4] or layer- 
ing of the SC material [5-7] . 

In this paper we consider a different route which may 
allow a substantial increase in T c for small SC wires, 
based on local cooling. To this aim we study a SC wire 
in contact with two different heat baths, held at different 
temperatures. We consider at first a wire in contact with 
two heat baths at its edges, held at different temperatures 
T L (left bath) and Tr (right bath), with Tr > T c e « > T L , 
where T^ q is the equilibrium critical temperature (up- 
per panel of Fig. 1(a)). We find that depending on the 
value of Tr,, the critical value Tr jC (defined as the max- 
imum value of the temperature Tr of the hot bath at 
which the wire is still SC) can be much larger than T^ q , 
which implies that the local temperature in the wire is 
maintained below T^ q although the average of Tr, and 
Tr exceeds it. We then study the effect of different cou- 
plings to the baths and of disorder on the above result. 
We find the dependence of Tr :C on the coupling to the 
different baths, and demonstrate that, in agreement with 
Anderson's theorem [8], weak disorder does not change 
Tr, c by much. Strong disorder leads to the breakdown of 
the SC state near the right (hot) bath, and to a spatial 
dependence of the SC order parameter. 

Normal metallic wires in contact with two heat baths 
at the edges were recently studied in detail, in the context 
of heat flow in such systems [9-11] . One of the main con- 
clusions of Ref. [10] was that in a system which is not dif- 
fusive, one cannot define a non-equilibrium temperature 
as the average of the temperatures of the different baths. 
Rather, it is the energy distribution function (DF) of the 



two baths which are averaged (a similar observation was 
verified experimentally for short wires, when interactions 
are relatively unimportant [12], a situation which is also 
likely to hold for the quasi-particles in SC wires). This 
observation will allow us to provide an analytic expres- 
sion for the critical temperature which shows excellent 
fit with our numerical calculations, and provides a direct 
prediction which may be tested experimentally. 

The method we use is a generalization of the 
Bogoliubov-De Gennes (BdG) mean-field theory [13] 
to non-equilibrium. The starting point is the tight- 
binding BdG Hamiltonian on a square lattice (with 
lattice constant o=l), Hsdc = E i0 -( e * ~ M) C L C «' ~ 
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ates an electron in the z-th lattice site with spin a, t 
is the hopping integral (t—1 serves as the energy scale 
hereafter), ej are random on-site energies drawn from 
a uniform distribution U[—W/2,W/2] (hence W is the 
strength of disorder, with W = representing a clean 
system), fi is the chemical potential and Aj is the SC or- 
der parameter in the z-th lattice site. The order parame- 
ter is to be determined self-consistently on every site via 
Aj = — £/(cjjCi|), where U > is the effective electron- 
electron attractive interaction and (•) stands for a statis- 
tical average. In this paper we treat only s-wave super- 
conductors, but the formalism can easily be extended to 
account for other kinds of symmetry. 

Since the BdG Hamiltonian is quadratic, it can be 
exactly diagonalized to describe the quasi-particle ex- 
citations j ntT . To diagonalize it, one performs a Bo- 
goliubov transformation [13] for the electron operators, 

4x = En ( u n(ihL + vKtfhnv), where u n (i) and v n (i) 
are quasi-particle and quasi-hole wave-functions, respec- 
tively. Requiring that this transformation diagonalizcs 
the Hamiltonian yields a set of eigenvalue equations for 
the quasiparticle (QP) wave functions [13], 
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where £u„(i) = -*X)<i,j,> + ( e * ~ A* + Uni/2)u n (i) 
(here nj is the local electron density), and Au n (i) = 
AiU n (i). From Eq. (1) one obtains the QP wave func- 
tions u n (i) and v n (i) and the energies E n . In equilibrium, 
this procedure results in a closed self-consistent set of 
equations for the local SC order parameter and density, 
which were recently used to study, e.g., effects of disorder 
and magnetic fields in two dimensional superconductors 
[14-16]. 

Since the QP excitations are non- interacting Fermions, 
they can be treated by the formalism of Refs. [9-11], 
which is aimed at studying such particles out of equilib- 
rium. To generalize the method of Refs. [9-11] to the 
QP excitations, we define a single-particle density ma- 
trix p, with matrix elements p nn > — {l n ^ln>-\) (note that 
the particle-hole symmetry allows one to treat only the 
up-spin excitations). The master equation for p is of 
the Lindblad form [17] (setting h — ks — 1 hereafter) 
p = —i{H,p] + £l[p] + £r[p] , where TL is the diagonal 
matrix of energies E n and 

4l,r)M = e K {e R)t e R \p} + e%e' R)t ) 

nn' 

(2) 

describe environment-induced inelastic transitions be- 
tween different single-particle states. The V-operators 
in Eq. (2) take on a local form [10] 



(L,R) 



,(L,R) /.(L,R) 



-,(L,R) 



r 



{u n {n)u* n ,{n) + v n (n)v*,(ri)) 



(3) 



where Sl,r are the contact area of the left 
(right) heat bath with the sample, /^' R ^(.E„) = 
1/(1 + exp(S„/TL j R)) is the Fermi distribution and To 
is some constant scattering rate. (We take T =0.1, 
changing r does not alter the results presented above.) 
Note that the V-operators operate on the QPs and 
not on the electron operators, since the Fermi function 
is defined for the occupation of the QPs (to put it 
differently, the part of the electrons which is in the 
superfluid phase does not feel scattering from the baths, 
only the QPs do). 

Once the V-operators are evaluated, the master equa- 
tion is solved in the asymptotic time limit (i.e., p = 0) 
and a solution for p is obtained. From this solution, the 
local density and order parameter are evaluated via the 
self-consistency condition, which reads 

A, = U^u n (i)v n (i)(l-2p nn ) 

n 

n t = 2j2(\ U n(i)\ 2 Pnn + \v n (i)\ 2 (l- Run)) ■ (4) 

n 

With Aj and nj determined, the whole procedure (i.e., 
the solution of the BdG equations, the evaluation of the 



V-operators and the solution of the master equation) is 
repeated until Aj and n-i no longer change (within the 
numerical tolerance of < 10 -5 ). We point that this treat- 
ment is of a mean-field type, and as such neglects statis- 
tical fluctuations in the Hamiltonian [18] or phase-slips 
[19]. As we will show below, in the ideal case one can 
define an effective temperature of the wire, and hence in 
such a case one can follow the usual treatment of phase- 
slips in SC wires, but using the effective temperature 
as input. Phase fluctuations are unlikely to change the 
local effective temperature, since these are excitations 
that do not carry heat (as opposed to QP excitations). 
In the disordered case, such a treatment cannot work as 
the temperature becomes space-dependent and one needs 
a different formalism to treat phase-fluctuations in the 
presence of a temperature gradient. Such a theory is be- 
yond the scope of the present work and will be the subject 
of future studies. 

We begin by presenting the averaged order parameter 
A = J2i A» , where V is the total number of lattice 
sites. We consider the geometry shown in the side panel 
of Fig. 1(a), where the SC wire (gray area) is connected 
to the thermal baths only at its edges. The numerical 
parameters are as follows. The wire dimensions are 100 x 
10, U = 2, W = (clean system) and the density is held 
at n = 0.875 (i.e., n electrons per site on average, and the 
chemical potential is chosen self-consistently to maintain 
this filling). We define Tj, by its ratio with the critical 
temperature at equilibrium T£ q , T^/Tj? 9 = 7 < 1. 

If Fig. 1 (a) , A is plotted as a function of Tr for differ- 
ent values of 7, 7 = 0.05, (top curve), 0.1, 1 (bottom 
curve). For low values of 7 , Tr jC can greatly exceed Tjf 9 
(marked by a solid arrow). In the inset of Fig. 1(a) we 
plot the local order parameter as a function of position 
along the wire (averaged over the transverse direction), 
at 7 = 0.05 for different values of Tr = 0.035, 0.385, 0.735 
and 3.185 (in units of t). We find that although there are 
two different temperatures at the edges, the order param- 
eter is practically uniform along the wire, in agreement 
with the results of Refs. [10, 11]. 

In Fig. 1(b) we study a somewhat different (and per- 
haps more realistic) situation, in which the heat baths 
are in contact with the wire not only at the edges but 
over some area (see side panel of Fig. 1(b)). We define 
the parameter a (1 — a) to be the ratio between the con- 
tact area of the right (left) heat bath and the area of the 
whole wire, such that a = 1 stands for a system in full 
contact only with the right heat bath. In Fig. 1(b), A 
is plotted as a function of Tr (at 7 = 0.083) for differ- 
ent values of a, a = 0.1, (top curve), 0.2, ...,1 (bottom 
curve). The a = 1 curve is the equilibrium curve (with 
T^ q marked by an arrow). As seen, for different values 
of a, Tr. c may again exceed T^ q . Also in this case we 
found that both A and the local temperature (calculated 
for a similar one-dimensional geometry, with the method 
of Refs. [fO, 11]) are uniform in space (not shown). This 
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emphasizes the fact that the temperature is defined not 
by the local baths, but rather by the (inelastic) scattering 
between states, which in the clean case span the entire 
system. 



T L /T*i = 7 we find 




FIG. 1: (color online) (a) Main panel: A as a function of 
Tr for different values of 7. At small values of 7, Tr can 
greatly exceed the equilibrium critical temperature (indicated 
by a solid arrow) without destroying superconductivity in the 
wire. Inset: position dependence of the order parameter for 
different values of Tr (at 7 = 0.05), demonstrating that it 
is uniform along the sample, (b) Same as in (a) but for a 
constant 7 = 0.083 and different values of a (describing the 
contact area of the right heat bath with the sample). Upper 
panel: the geometries considered in (a) and (b) . 

If one could simply define a local temperature which 
gradually shifts from Tl to Tr, then one would expect 
that Tr could not exceed T^ q and that the order pa- 
rameter would not be uniform in space. In order to ex- 
plain our findings, we recall that one of the main results 
of Refs. [10, 11] is that in the ballistic limit the tem- 
perature is uniform, and a non-equilibrium DF develops, 
which is the average of the two DFs of the left and right 
baths. In the case represented in Fig. 1(b), we find that 
a weighted average between the DFs of the left and right 
baths develops, the weight being a. This result, along 
with the observation that Aj is uniform in space, allows 
us to find an analytical expression for the effective T c as 
follows. In the equilibrium theory of superconductivity 
[13], the critical temperature T^ q is determined by the 
gap equation = \ ^ D,Tc " i(l - 2f(x))dx ,where 
N a is the density of states at the Fermi energy, ujy 
is the Debye frequency and f(x) is the DF. From the 
above discussion, in the non-equilibrium case we have 
f(x) — afjj(x) + (1 — a)fp(x). The resulting equa- 
tion for Tr c then reads 



N U 



-(1 - 2f(x))dx + 
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These integrals may be evaluated exactly, and with 



Tr.c 



(6) 



In Fig. 2 we plot T R>C /T^ as a function of 7 (Fig. 2(a)) 
and of a (Fig. 2(b)), taken from the data of Fig. 1(a) 
and(b), respectively. The solid line corresponds to 
Eq. (6) for the two cases, with the corresponding param- 
eters taken from the numerical calculation. The agree- 
ment between Eq. (6) and the numerical results confirms 
that indeed the DF is a weighted average of the DFs of 
the two baths. One can now use Eq. (6) to estimate the 
effective Tr^ c of real materials. As a practical use, it is 
advantageous to raise Tr iC above the freezing point of 
liquid Nitrogen, ~ 77K . For example, consider a desired 
working temperature of ~ 80K, and local refrigerators 
which cool down to 40K , deposited on a SC wire with 
Tj? 9 ~ 60K. A cover of 40% refrigerators would increase 
T fl , c to = 78K. 

Perhaps a more intriguing possibility is the enhance- 
ment of Tr c to room temperature. For a wire made 
of the newly-found Iron compound [20-22] (T c e « - 50K) 
heated (or cooled) at the edges (a = 0.5), a temperature 
Ti = 7.5K would drive Tr c above room temperature. 
For nano-scale wires made of a high-T c material, the fab- 
rication of which was recently demonstrated [23], taking 
T^ ~ 80K, local cooling of T L = 20K and coverage of 
a = 0.5 would drive Tr, c above room temperature. Of 
course, other effects (e.g. phonon scattering, phase fluc- 
tuations etc.) might become very important in such high 
value of Tr and inhibit the zero-resistance state. 
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FIG. 2: (a) The ratio TR, C /T" q as a function of 7, taken from 
the data of Fig. 1(a). The points correspond to the numerical 
data and the solid line is Eq. (6) with a = 0.5. (b) Tr iC /T^ 
as a function of a, taken from the data of Fig. 1(b). The 
points correspond to the numerical data and the solid line is 
Eq. (6) with 7 = 0.083. 

Next we turn to study the effect of disorder. In 
Rcf. [11] it was shown that the form of the non- 
equilibrium DF is robust against disorder, but that 
the local temperature profile changes from a constant- 
temperature to a position-dependent profile. In Fig. 3 
we plot the order parameter A^ as a function of position 
along the wire for Tr — 1.5 (which is above T^ q ~ 0.2 
in this example), Tl — 0.02 and dimensions 50 x 10, for 
different values of disorder, W = (dark curve), 0.2, 4 
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(bright curve). The order parameter is averaged over 500 
realizations of disorder. For a clean system A is uniform, 
and assumes a position-dependence with increasing dis- 
order. For small values of disorder it is finite everywhere 
in the sample, in agreement with Anderson's theorem 
[8] . For large values of disorder it vanishes near the right 
edge, and increases at the left edge. This is due to the fact 
that a local temperature ensues which varies from the left 
to the right temperature. In fact, if one compares Fig. 3 
to the local temperature calculated for similar (normal) 
systems (Fig. 1 of Ref. [11]) one finds that A obeys a 
simple BCS-like law, with the temperature replaced by 
the local temperature. 

We point out that in this case the local temperature 
near the left edge (which corresponds to Tl) is lower than 
the effective temperature of the clean sample, and hence 
the rise in the value of A with increasing disorder near 
the left edge. The distance from the right edge at which 
A vanishes indicates the "thermal length" where the lo- 
cal temperature is close to Tr [11]. We also note that 
already for relatively small values of disorder (W « 0.5 
in our case) the localization length is smaller than the 
system size, which means that the onset of a vanishing 
gap in the system (occurring at W w 1.4) does not di- 
rectly correspond to the onset of localization [11, 24]. 
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FIG. 3: Position dependence of the order parameter A, cal- 
culated with increasing values of disorder, from the clean case 
W = (dark line) to strongly disordered case W — 4 (light 
gray line), averaged over 500 realizations of disorder. For this 
example T L = 0.02, T^ q = 0.2 and T R = 1.5. Inset: Sug- 
gested experimental setup. An insulating layer is deposited 
on top of a SC wire with etched contacts, and a heater coil 
is placed on top of it, to generate local heating. The resis- 
tance is then monitored along the other contacts (which may 
be either 4-terminal contacts or tunnel junctions) . 

In order to test our predictions, we suggest the experi- 
mental setup shown in the inset of Fig. 3. It consists of a 
SC wire with etched contacts (either regular 4-terminal 
contacts, or made as tunnel junction, aimed at measuring 
the local density of states at that location). On top of 
one of the contacts an insulating layer is deposited, and 
on top of it a heater coil is set. The temperature below 
the heater can be calibrated by measuring the resistance 



between the contacts beneath the heater when the wire 
is in the normal state. Then the whole device is cooled 
down, and the resistance through the other contacts is 
measured as a function of the current that passes through 
the heater (i.e., the local temperature beneath it). In a 
uniform system, the resistances of all the contacts should 
vanish if the system is SC, and have finite values once the 
local temperature beneath the heater rises above Tr jC . If 
the system is disordered, the different contacts should 
exhibit a gradual transition to a normal state. 

The length-scale that determines the onset of a tem- 
perature gradient is in this case the QP mean-free path 
(or diffusion length) [24]. While it can be estimated for 
low T c metals, for high-T c materials it is unknown (al- 
though it is suspected to be small, based on their poor 
conduction in the normal state). By controlling the dis- 
tances between the contacts our proposed experiment can 
thus serve to determine this length, by relating it to the 
length-scale at which a temperature gradient develops. 

In this work we have neglected the effect of phonon 
scattering, and assumed the effective electron-electron 
interaction is not appreciably changed by temperature, 
supported by the fact that the Debye temperature is 
much larger than T c . However, considering our geometry, 
electron-phonon (e-ph) interaction effects may play a sig- 
nificant role as Tr reaches the Debye temperature. This 
is probably more important in disordered wires (for clean 
wires we expect that, on equal grounds, phonons will 
also acquire a uniform temperature). The e-ph interac- 
tion may induce inelastic electron transitions, which will 
reduce the inelastic mean-free path. Since the temper- 
ature profile is sensitive to the inelastic mean-free path 
[24], this effect may be seen in our suggested experiment. 

To conclude, we point out that in recent years there 
have been tremendous advances in fabricating micro- 
refrigerators, based on the thermo-electric Peltier effect, 
and which can locally cool down their environment sub- 
stantially [25, 26]. Since the efficiency of thermo-electric 
materials is likely to increase in the future [27], one can 
conceive a device composed of a (relatively) high-T c ma- 
terial, on top of which are embedded a series of micro- 
refrigerators, covering an area of the material and cool- 
ing it enough for it to operate at a temperature which 
exceeds its T c , a possibility that may allow for integrat- 
ing superconducting wires as circuit elements in various 
devices. 
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